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Context. The blazar 1ES 1101-232 was observed with the High Energy Stereoscopic System (H.E.S.S.) of Atmospheric Cherenkov 
Telescopes (ACT) in 2004 and 2005, for a live time of 43 hours. VHE (£>10 n eV) 7-rays were detected for the first time from 
this object. 

Aims. VHE observations of blazars are used to investigate the inner parts of the blazar jets, and also to study the extragalactic 
background light (EBL) in the near-infrared band. 
• i-h . Methods. Observations in 2005 were conducted in a multiwavelength campaign, together with the RXTE satellite and optical 
' observations. In 2004, simultaneous observations with XMM-Newton were obtained. 
$H ' Results. 1ES 1101-232 was detected with H.E.S.S. with an excess of 649 photons, at a significance of 10 a. The measured VHE 
5^ ! 7-ray flux amounts to dN/dE = (5.63 ± 0.89) x 10 _13 (-B/TeV) _(2 94±0 ' 20) cm~ 2 s~ 1 TeV~ 1 , above a spectral energy threshold of 
225 GeV. No significant variation of the VHE 7-ray flux on any time scale was found. 1ES 1101-232 exhibits a very hard spectrum, 
and at a redshift of z — 0.186, is the blazar with the highest confirmed redshift detected in VHE 7-rays so far. 
Conclusions. The data allow the construction of truly simultaneous spectral energy distributions of the source, from the optical to 
the VHE band. Using an EBL model with vF u = 14nWm~ 2 sr~ 1 at 1.5 pm as presented in Aharonian et al. (2006a) suggests an 
intrinsic VHE power output peak of the source at above 3TeV. 
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1. Introduction 

Blazars (BL Lacs and Flat Spectrum Radio Quasars) are 
thought to be active galactic nuclei (AGN) that have their 
jet axis oriented close to the line of sight of the observer. 
The broadband spectral energy distribution (SED, vF v 
representation) of blazars is characterized by two peaks, 
one at optical to X-ray energies, and another at 7-ray 
energies. The low-energy branch is commonly explained 
as electron synchrotron emission. The high-energy branch 
can be explained in a variety of ways. In leptonic scenar- 
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ios, it is assumed to result from Inverse Compton (IC) 
emission from the same electron population, up-scattering 
the self-generated synchrotron photons or external pho- 
tons (synchrotron self-Compton, SSC, and external 



vations are reported. In Section we concentrate on the 
interpretation of the spectral energy distribution derived 
from the source. 
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ternative hadronic scenarios, accelerated protons are 
the main source of high-energy radiation, either di- 
rectly or through the producti o n of s econdary particles 
(e.g.. | Biermann fe Strittmatterl 1987 ; Aharonian] 2000; 



iPohl fe Schlickeiserl l2000t iMiicke fc Protherod l2001ft . 
Observationally, the SED and the variability in the dif- 
ferent bands carry the information about the acceleration 
processes at work in the jet, and could ultimately also 
shed light on the energy transfer mechanism of the central 
engine - a supermassive black hole - into the jet. 

Over the past fourteen years, VHE 7-ray emission from 
appr oximately a dozen blazars has been detected (see, 
e.g., lOnd 120051 for a recent review). Both the detection 
of fast variability and the availability of broadband ob- 
servations - especially including X-ray measurements - 
have bee n used to constrain individual source parame - 
ters (e.g.. lKrawczynski et al.ll2002l ; I Aharonian et al.l2005aft . 
Increasing the number of known VHE blazars, especially at 
higher redshift, is of importance for two reasons: 

(1) Relatively little is still known about the average be- 
haviour of VHE blazars. Most VHE blazars detected so 
far belong to the classes of X-ray selected BL Lacs (XBL) 
or high-frequency peaked BL Lacs (HBL), but popula- 
tion studies are restricted by the low number of sources. 
Previous detections of VHE blazars have also been biased 
towards high states of the sources, because of the lim- 
ited sensitivity of the available instruments. It was shown 
only recently that quiescent sta tes can be detected now 
in short (~hours) observations (jAharonian et al.ll2005aft . 
Little is known about average activity cycles and flare 
time scales, except for a few sources: Mkn 421 (e.g., 
lAharonian et alj l2002bt iBlazeiowski et al.ll2005ft. Mk n 501 
(e.g.. lAharonian et all 1999: Krawczynski et al. 2002) , and 
1ES 1 959+650 (e.g.. iKrawczvnski et alj|2004 lAlbert et all 
12006ft . 

(2) Source photons above ~100 GeV are attenuated 
by the EBL through 7-7-interactions. Therefore, an EBL 
density in the relevant waveband range (typically ~0.1 
to ~10/im) has to be assumed to derive the intrinsic 
blazar spectrum. Conversely, if it is possible to determine 
or constrain the intrinsic blazar spectrum through mod- 
els, then constraints on the EB L density in the respective 
waveband can be obtai ned (e.g. ICoppi fc Aharonian] 1 19991 : 
lAharonian et al.l [2006a] hereafter AHA06a). 

To date, 1ES 1101-232 is the most distant VHE blazar 
known with confirmed redshift (z — 0.186). It should be 
noted that the hard spectrum observed from 1ES 1101-232 
and its relatively large redshift allow for strong constraints 
on the EBL dens ity. This is described in detail in a separate 
paper (|AHA06al ). 

In this paper, the discovery of VHE 7-ray emission from 
the blazar 1ES 1101-232 with H.E.S.S. is reported. The pa- 
per is organized as follows: In Section [2] the source char- 
acteristics of 1ES 1101-232 are described. The results of 
the H.E.S.S. observations of 1ES 1101-232 are presented in 
Section [3] In Section [4l multifrequency observations that 
were performed contemporaneously to the H.E.S.S. obser- 



Emission from 1ES 1101-232 was first detected by the Ariel- 
5 X-ray satellite, the source A 1059-22 was however misiden- 
tified with th e Abell 1146 cluster of galaxies at z=0.13 9 
at that time (|McHardv et alJll98lt iMaccagni et al.lll978l ). 
The HEAO-1 source H 1101-232 was later correctly identi- 
fied as a BL Lac type object, using the optical and radi o 
counterparts ([Buckley et all 119851 : IRemillard et"all 11989ft . 
The source has been detected by multiple X-ray observa- 
tories, and for the purpose of this paper, the commonly 
used na me 1ES 1101-232 from the Einstein sle w survey is 
adopted ([Elvis et alJll992UPerlman et al.lll996ft . 

1ES 1101-232 resid es in an elliptical host galaxy at a 
redsh ift of z = 0.186 (|Remillard et all 119891 : IFalomo et all 
11994ft. The host is presumably part o f a galaxy cluster 
(IRemillard etaHl 19891: IPesce et al.lll994ft . VLA maps of the 
BL Lac show a one-sided diffuse structure to the north of 
<~ 45" size, but no wel l-collimated jet outside a few kpc dis- 
tance from the core (|Laurent-Muehleisen et aLlll993ft. In 
the optical, the host galaxy is resolved (Rcmil lard et all 
119891 : lAbraham et al.lfl991ft . IFalomo et aTJ (|1993ft deduced 
its brightness using a spectroscopic deconvolution of host 
and BL Lac. The most recent esti mate of wr = 16.41 wa s 
derived from an ang ular profile fit (jFalomo fc Ulrichl feOOO) . 
The galaxy is one of the brightest BL Lac host galaxies 
so far detected, and also the only one known with sig- 
nificantly boxy isophotes (|Falomo fc Ulrichl feOOO). indicat- 
ing a merger process or extra dust components. The BL 
Lac itself has a typi cal brightness of my — 16-17 (e.g., 
IRemillard etall 1989ft . The optical emission from 1ES 1101- 
232 has typically var ied on the timescale of months (e.g., 
IRemillard et al .Ml 989ft . Optical flares on intra day timescales 
have also been claimed in one observation (|Romero et all 
11999ft . 

The source has been classified earlier as an XBL (e.g., 
Scarpa fc Falomol 1997ft . and later on as an HBL (e.g. 



Donato et al. I 1200 lft~ because of the dominance of syn- 



chrotron emission in the X-ray band. Several authors have 
concluded from the broadband characteristics of 1ES 1101- 
232 that this source is expected to emit VHE 7-ray emis- 
sion at flux levels detect a ble by instruments like H.E.S.S . 
(e.g., IWolter et all l2000t ICostamante fc GhisellinJ 12002ft . 
Previous VHE observations with the Durh am Mark 6 tele- 
scope in 1998 have only yielded flux limits ([Chadwick et al.l 
1999). Also, in the GeV 7-ray do main, EGRET d id not de- 
tect emission from 1ES 1101-232 (|Lin et aLjll996ft. 

In previously published SSC models (jWolter et alj|2000t 
ICostamante fc Ghisellinil 12002) . the IC peak was generally 
expected to be around 100 GeV, but this seems not to be 
the case as shown in this paper. 



3. H.E.S.S. observations of 1ES 1101-232 

Observations were made with the H.E.S.S. Cherenkov tele- 
scopes in April and June 2004, and in March 2005. On June 
8 th 2004, also XMM-Newton X-ray observations were ob- 
tained, scheduled such that simultaneous H.E.S.S. observa- 
tions could be conducted. Following the detection of a weak 
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signal in the 2004 H.E.S.S. observations, an extended mul- 
tifrequency campaign was organized for 11 nights in March 
2005, to study the broadband emission from 1ES 1101-232 
and to search for (possibly correlated) variability in the 
different wavebands. Simultaneous observations were car- 
ried out with H.E.S.S., X-ray measurements with RXTE, 
and optical measurements with the ROTSE 3c robotic tele- 
scope. 

3.1. The H.E.S.S. experiment 

The H.E.S.S. experiment (jHinton et al.l l2004h consists in 
phase I of four ACTs operating stereoscopically. Each tele- 
scope consists of a tesselated 13 m diameter (107 m 2 surface 
area) mirror, which focuses the Cherenkov light from the 
showers of secondary particles created by the interaction of 
7-rays in the atmosphere onto a camera in the focal plane. 
Each camera consists of 960 photomultipliers with a pixel 
size of 0.16°, providing a field of view (FoV) of 5°. The array 
is located in the Khomas highlands in Namibia (—23° 16', 
16°30', 1835m a.s.L). 

The angular resolution of the stereo system is better 
than 0.1° per event. The energy threshold of H.E.S.S. is 
about 100 GeV (at zenith), with spectral measurements 
possible above ^150 GeV with an energy resolution of 15%. 
The energy threshold increases with zenith angle. For the 
data set of 1ES 1101-232 discussed in this paper, the time- 
averaged spectrum presented in this paper has an energy 
threshold of 225 GeV. 

The minimum detectable point source flux above 
100 GeV with H.E.S.S. is ~ 4 x lO^ei-gcm^s" 1 for a 
5cr detection in 25 hou rs, corresponding to ^TOmCrab 
([Aharonian et al.ll2006cf ). The sensitivity enabled a > 3 c 
detection per night (~ 5 hours) in the case of the 1ES 1101- 
232 observations. 
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Table 1. Log of the H.E.S.S. observations on 1ES1101- 
232 in 2004 and 2005. Numbers reflect the observa- 
tions after data cleaning and good run selection. Nights 
with observations performed on 1ES 1101-232, where how- 
ever all H.E.S.S. observations needed to be discarded 
because of weather selections, are marked with a '- 
'. 1 Modified Julian date. 2 Live time T[hours]. 3 Flux 
F(£;>200GeV)[10- 12 ergcm" 2 s~ 1 ]. 4 Detection significance 
S in units of standard deviations. 5 Upper limit at 99% con- 
fidence level. 



3.2. H.E.S.S. observations of 1ES 1101-232 

First observations of 1ES 1101-232 with H.E.S.S. were per- 
formed during four nights in April 2004, for a total live 
time of 2.7 hours after quality selection, and for six nights 
in June 2004, for a total of 8.4 hours after quality selection. 
The total data set in March 2005 after quality selection 
amounts to 31.6 hours live time. The observation log is 
shown in Tab. [I] 

H.E.S.S. observations were taken in runs of typically 
28 min each. Runs were taken in wobble mod^B, with a dec- 
lination or right ascension offset of ±0.5° in 2004 and ±0.7° 
in 2005. Since the 7-ray acceptance of the H.E.S.S. instru- 
ment is nearly constant within 1° radius from the FoV cen- 
ter, this change of observation mode results in a slight sen- 
sitivity increase, as a larger background control area with 
equal acceptance in the FoV can be used. 

The data cleaning to derive a set of good quality runs, 
which are used in the data analysis, consists of two filter- 
ing processes. First, from the individual shower images as 
recorded from triggered telescopes, all pixels that have not 
worked properly are removed; occasionally, entire cameras 
have to be excluded from individual runs. Then runs are 



1 In wobble mode, the source is displaced with respect to 
the center of the FoV, the sign of the offset alternating be- 
tween consecutive runs. This permitted continuous monitoring 
of 1ES 1101-232. 



discarded that show a too low or fluctuating event trigger 
rate, caused by bad atmospheric conditions or clouds. 
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Fig. 1. Angular event distribution. Events are from the en- 
tire H.E.S.S. data set on 1ES 1101-232, after shape cuts 
to suppress a large fraction of the background. The filled 
circles denote the event distribution in squared distance 
with respect to the direction to 1ES 1101-232. The filled 
histogram represents the expected background. For this 
histogram, the distributions obtained from seven control 
regions were summed and renormalized. 
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Table 2. H.E.S.S. analysis parameters and event numbers. 
Values are given for the total H.E.S.S. data set as well as 
for the three data subsets considered. 



3.3. Analysis of the H.E.S.S. data 

The H.E.S.S. data have been processed and analysed ac- 
cording to the standard an a lysis ch ain a s detailed for in - 
stance in lAharonian elaTI (|2005bJ ) and iBenbowl (|2005f) . 
Shower images as recorded by individual telescopes are 
parametrized using the Hillas parametrisation. The shower 
direction (i.e., the origin of the incoming 7-ray) and the pro- 
jected shower impact point on the ground are determined 
from a stereoscopic reconstruction of the air shower. A large 
fraction of the background is suppressed using cuts on the 
shower shape parameters mean reduced scaled width mrsw 
and length mrsl. Photon directions are used for a one- 
dimensional projection including background estimate (9 2 - 
plot, Fig.[T]) and to reconstruct the source location (Fig. [2]). 
The background after shape cuts is estimated from different 
control regions in the FoV, having the same distance to the 
center of the FoV i n order to be indep endent of the radial 
acceptance change (jHinton et al.ll2005T) . For this data set, a 
background region with an area A s 11 (in 2004) or 15 (in 
2005) times larger than the on-source area A on was used. 
The applied standard cuts (see Tab. [2]) were optimized on a 
simulated source with 10% of the flux from the Crab Nebula 
and a Crab-like power-law spectrum oc E~ r with T=2.6. 

Energy estimates for individual photons are based on 
the comparison of the measured image amplitudes to data 
from simulated events, using the measured shower core po- 
sition and zenith angle as parameters. The resolution per 
event is ~15% above the spectral threshold, which is de- 



fined as the energy threshold above which the energy re- 
construction bias is less than 10%. Only events above this 
safe threshold are used to compute energy spectra. The ef- 
fective area is based on Monte Carlo simulations and is 
depending on zenith angle and system configuration (i.e., 
which telescopes are included in a particular run). 

Compared to the analysis results used in AHA06a, 
an improved energy calibration of the telescope system 
was applied to the data, better taking into account the 
long-term optical sensi tivity changes of the instrument 
(Ah aroman et al.ll2006cl ). For the given total data sample, 
this energy scale recalibration yields a safe energy th reshold 
of 225 GeV (compared to 165 GeV used in lAHA06aD and a 
flux normalisation increase of 27% at 1 TeV. After this cor- 
rectio n, the systematic flux un certainty is now estimated as 
20% (| Aharonian et al.l l2006d ). Reconstructed spectral in- 
dices were not affected significantly by these calibration 
updates, the systematic error estimate for reconstructed 
photon indices is Ar sys ~ 0.1 ( AHA06ai; lAharonian et all 
(20063)). The recalibration slightly increased the back- 
ground noise in this data set which therefore has an excess 
significance of 10.1 cr, slightly smaller than the detection sig- 
nificance of 11.7 cr th at was derived from the original data 
set used in lAHA06al . 

The results derived have been verified using a com- 
pletely independent calibrat ion and analysis chain, w hich 
is described for insta nce in lAharonian et al.l (|2006bl ) and 
iLemoine et "ail (|2006[ ). 

3.4. Skymap and identification of the VHE "/-ray source 

Source confusion is generally a minor issue in VHE data, 
especially from extragalactic sources. No other source 
near 1ES 1101-232 (ll h 03 m 37 57 s , -23°29'30.2'', J2000, 
lLaurent-Muehleisen et al.l I1993D is a good candidate for 
VHE 7-ray emission. The VHE 7-ray source location 
was derived from the skymap of VHE photons and is 
ll h 03 m 36.5 s ± 2.5 s tat , -23°29'45" ± 26£ at (J2000), which 
is consistent with the above given radio position, and also 
with the optical and X-ray positions as shown in Fig. [21 As 
the present VHE data do not exhibit variability, the iden- 
tification of the VHE 7-ray source is based on its location 
and the interpretation of the SED. 

Figure [2] also shows that the X-ray imaging data, which 
were obtained with XMM-Newton, are compatible with the 
XMM-Newton point spread function. This excludes the in- 
fluence of possible nearby X-ray sources in the FoV of non- 
imaging X-ray measurements such as RXTE. 

3.5. VHE j-ray light curve 

The fluxes derived from the three H.E.S.S. data sets (April 
2004, June 2004, March 2005) are statistically compati- 
ble with each other, indicating that the VHE 7-ray flux 
has remained constant throughout these observing periods. 
However, as the detections from the April (1.9 a) and June 
2004 (3.7 cr) data sets alone are not very significant, only 
flux variations of a factor of approximately two or larger can 
be excluded. Night-by-night variations were also searched 
for, but the corresponding light curve is compatible with a 
constant flux (probability for a constant flux of the total 
data set 85%, of the March 2005 data set 64%). 
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1 l h 03 m 40 s 



1 l h 03 m 32 s 



R.A. (2000) 



Fig. 2. Sky map of the region surrounding 1ES 1101-232. 
An R-band image made by the Anglo-Australian obser- 
vatory with the UK Schmidt telescope is shown in grey- 
scale. The host galaxy (labeled 1ES 1101-232 in the im- 
age) of the BL Lac is resolve d as an elliptical galaxy , with 
boxy isophotes at larger radii (|Falomo fc UlricMl200a . The 
dashed ellipse denotes the 1 a error of the reconstructed 
VHE 7-ray source position. The solid contours are from the 
XMM-Newton X-ray measurements with the MOS 2 cam- 
era, as discussed in Section [4721 The contour spacing is log- 
arithmic (0.1, 0.3, 1.4, 5.8 and 24% of the peak intensity), 
the outermost contour also denotes the edge of the effective 
MOS 2 FoV for this small window mode observation. 



In Fig. [3l the nightly averaged light curve of H.E.S.S. 
is shown for the March 2005 period, together with si- 
multaneous X-ray measurements from the RXTE satellite. 
The VHE 7-ray flux is given in integral flux units above 
200 GeV and was computed under the assumption of the 
time-averaged spectral index of T — 2.94. The X-ray fluxes 
were similarly derived, details are given in Section 14.1.11 
The RXTE data indicate an X-ray flux variation of ~15% 
(min-max), whereas the simultaneously taken H.E.S.S. data 
are not sensitive enough to detect possible correlations with 
similar amplitudes in the X-ray and VHE bands. 

3.6. VHE j-ray spectrum 

The differential energy spectra derived for the entire 
H.E.S.S. data set, as well as for the June 2004 and the 
March 2005 data subsets, are shown in the left pan- 
els of Fig.[4l The low statistics of the April 2004 data 
set prevented us from performing a spectral analysis on 
that subset. The measured spectra are compatible with 
power laws, Tab. [3] lists the corresponding photon indices 
r a bs derived from fits between 0.2 and 4TeV. The fit 
for the total spectrum gives dN/dE = (5.63 ± 0.89) • 
10- 13 (£;/TcV)-( 2 - 94±0 - 20 )cm- 2 s- 1 TeV" 1 above a spectral 
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Fig. 3. March 2005 light curve of 1ES 1101-232. 
MJD=53435 corresponds to the night Mar. 5.-6. Upper 
panel: VHE 7-ray flux as measured with H.E.S.S., in 
nightly averaged bins. Within errors, the flux was con- 
stant. Lower panel: X-ray flux (2-10 keV), derived from 
the PCU2 detector onboard RXTE. Note that the flux 
scale is truncated, the X-ray flux difference from minimum 
to maximum is ^15%. The thick horizontal bars in both 
panels denote the times (first to last measurement of an 
observing night) when the VHE and X-ray data were 
taken, indicating the high degree of simultaneity of the 
two data sets. 



energy threshold of 225 GeV. The integral flux is F(E > 
200 GeV) = (4.5 ± 1.2) ■ lO-^ergcm-^" 1 . 

Spectral bins above 4TeV do not contain a significant 
flux detection. Therefore, from these bins 99% upper lim- 
its were derived and shown in Fig.[3J The photon indices 
Tabs were derived excluding these bins. Table [3] addition- 
ally lists photon indices r* bs from single power law fits, 
for which these flux estimates above 4TeV were included 
in x 2 -fits. This is a viable procedure since the error bars 
are dominated by background estimates and are therefore 
mostly Gaussian. However, the results depend on the as- 
sumption that the power laws would extend to these high 
energies. 
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3. 7. VHE spectrum corrected for EBL absorption 

In the following, the data sets were corrected for EBL 
absorption, to investigate whether the intrinsic spectra 
show evidence for a cut-off towards high energies. Spectral 
changes between periods were also searched for. The deab- 
sorbed spectra are used in the SED, as discussed in Section 

El 

The intrinsic spectrum of 1ES 1101-232 depends on the 
assumed EBL spectrum. For a detailed discussion about 
the EBL abso rption of VHE 7 -ray spectra, we refer the 
reader to, e.g.. lAharonianl (|2001h . Corrections with a range 
of plausible EBL spectra result in deabsorbed spectra of 
1ES 1101-232 that can be described over the entire detected 
energy range, i.e. between 0.2 and 4TeV, by a single power 
law of photon index r deab s . i.e., r deabs = r abs - AF, see 
lAHA06at IStecker fc Scullvl (|2006tl . 

3.7.1. Maximum EBL 



To represent what AHA06a considered as th e highest pos- 
sible EBL level, the EBL shape P0.45 from |AHA06a| was 
used (cf. also Section [3~8f . This shape reflects an EBL level 



of 14nWm _2 sr _1 at 1.5 /mi after scaling down by 15% to 
take galaxy evolution effects into account. The such de- 
absorbed spectra of 1ES 1101-232 are shown in the right 
panels of Fig.UJ The fit to the deabsorbed spectrum of 
the total data set in the energy range 0.2 to 4TeV yields 
Tdcabs = 1.51 ± 0.17. The spectra of the two subsets consid- 
ered (r Jun200 4 = 1.70 ± 0.47 and r M ar2005 = 1-49 ± 0.19) 
are statistically compatible with each other. 

Similar results are obtained when using the EBL 



baseline" model described in Stecker et al.l (|2006f ). as 



sh own bv IStecker fc Scullvl (|2006|). Also the EBL model 



by iPrimack et al.l ( 20011 ) yields similar numbers, after a 
slight reduction of 15% to match the above quoted level 
of 14nWm _2 sr _1 . 

It is of interest to test whether the upper limits above 
4 TeV are compatible with a power-law extrapolation of the 
lower energy spectra, or are indicative of a steepening of the 
spectrum. For the intrinsic spectra, this would imply an ob- 
servational hint for a peak in the SED at that energy range. 
For the deabsorbed spectra, only those bins above 4 TeV 
could be included where the chosen EBL parametrisation 
(P0.45) can still be considered valid, which restricts the 
range to < 7 TeV. As justified in Section [331 x 2 -tests were 
applied. For example, for the March 2005 spectrum, the 
probability of the spectrum being compatible with T = 1.49 
changes from 65% to 49%. Hence, there is no observational 
hint for a steepening of the spectrum above ~ 4 TeV, but 
the data do not exclude a spectral change above that en- 
ergy either. Photon indices r*j cabs from power-law fits in 
the energy range 0.2 to 7 TeV are listed in Tab. [31 



3.7.2. Further EBL levels 

Lowering the EBL level used for the deabsorbtion 
leads to a softening of the spectra. Lower limits in 
the rel evant EBL waveband r ange come from galaxy 
counts (iMadau fc PozzettH l2000t ) and are of the order of 
10nWm~ 2 sr~ 1 at 1.5 jtim. The resolved galaxy counts may 
however represe nt a too low EBL level because of possi- 
bly missed light (jTotani et al.ll200i[ ). Nevertheless, to rep- 
resent the such constrained minimum EBL, the representa- 



tion P0.40 ((Ah aroman "et"a!ll2006a|) was chosen and scaled 
down by 15% to take galaxy evolution effects into ac- 
count. Using this minimum EBL to deabsorb the spectra 
result in Tdoabs = 1-85 ± 0.18, Fj un2 oo4 = 2.05 ± 0.56 and 

r M ar2005 = 1.84 ± 0.20. 

EBL models hig her than the maximu m EBL level were 

" (|2006l . 



described, e.e 



Stecker et al.l 



"fast evolution" 



case) and by frCneiske et al. I (|2004 "best fit" model), with 
an EBL den sity of about uF y ( 1.5 a m) ~ 20nWm~ 2 sr 



As shown in IStecker fc Scullvl (j2006t ), the "fast evolution" 
EBL would lead to an intrinsic spectrum with Tdoabs — 1-0. 
This result would be in conflict with the assumption of a 
limit on the intrinsic hardness of VHE blazar spectra, see 
next Section and l5.ll 

3.8. EBL limit 



Following the discussion in lAHA06al we assume that the 
intrinsic blazar spectrum did not have a photon index sig- 
nificantly harder than 1.5 (i.e., Tdcabs < 1-5), taking the 
present observational and theoretical knowledge of VHE 
blazar spectra into account. Applied to the data from 
1ES 1101-232, this limit results in an upper limit of the 
EBL density in the waveband range of 0.75-3.5 /mi, with a 
peak density of (14 ± 4) n Wm~ 2 sr~ 1 at 1.5 /mi for typical 
EBL shapes as reported in lAHA0"6al 

Given that the updated calibration of the data set 
(Section I3.3[) involves a change of the energy calibration 
which is s lightly lar ger than the systematic error of ±15% 
quoted in AHA06a, the procedure described in that paper 
to derive the EBL upper limit was repeated. The updated 
analysis constrains the total spectrum to a power law with 
T = 2.94 between 0.23 and 4.0 TeV, whereas in lAHA06al 
T = 2.88 between 0.16 and 3.3 TeV was used. The EBL 
limit de rived with these updated numbers differs from the 
value in lAHA06a[ by less than 1%, well below any statis- 
tical uncertainties involved. The upper limit for the EBL 
spectrum as given above is therefore confirmed. 

The error quoted for the peak density mostly comes 
from the statistical error of the photon index derived 
from the 1ES 1101-232 data set. This uncertainty is there- 
fore well represented by the statistical errors of the deab- 
sorbed spectra of 1ES 1101-232 discussed in the previous 
Section. It is therefore inappropriate to translate the value 
of (14±4) nWm^sr" 1 into an 68% or 95% upper limit. We 
remark that - because of the procedure described above - 
the choice of 14nWm~ 2 sr -1 yields by construction an in- 
trinsic photon index of the total 1ES 1101-232 data set of 
1.5. 



4. Broadband observations of 1ES 1101-232 

4.1. The observation campaign in March 2005 
4.1.1. X-ray data 

110.2ksec RXTE observation time for monitoring obser- 
vations have been obtained, which were scheduled such 
that simultaneous H.E.S.S. observations were possible for 
11 consecutive nights in March 2005 (see Tab.[l|. Because of 
Earth occultation and downtime in the southern Atlantic 
anomaly (SAA), the average on-time of RXTE was 56% 
during the "nightly" satellite observation windows. 
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Fig. 4. VHE 7-ray spectra from 1ES 1101-232. Left panels: Reconstructed photon flux, as measured with H.E.S.S. 
Right panels: Photon spectra after correction for maximum EBL absorption, using an EBL model with 14nWm~ 2 sr _1 
at 1.5 (im as described in Section [3.7.11 Upper limits in these deabsorbed spectra at energies above 7 TeV are shown as 
open symbols only, because of strong EBL uncertainties at these high energies. In all panels, solid lines denote power-law 
fits between 0.2 and 4 TeV. Extrapolations of these power laws to higher energies are shown as dashed lines. 



RXTE data were analyzed using standard reduction 
routines. During most observations, two PCA detectors 
(PCU0 and PCU2) were active. For the analysis pre- 
sented here, only results from PCU2 were used. PCU0 
has lost its front veto layer, and is therefore suscep- 
tible to unmode led and unremoved background events 
(|Xue fc Cuill2005f ). Only the signal from the top layer (X1L, 
X1R) was used for optimum signal to noise ratio. The 
STANDARD2 data were extracted u sing XDF and processed 
using f tools (jBlackburnl 119951 ) from HEASDFT 6.0.4. 
The data were filtered using standard criteria recom- 
mended by the RXTE guest observer facility (avoidance 
of Earth's limb, SAA, and tracking offsets). Electron con- 
tamination was suppressed by limiting the corresponding 
ELECTRDN2 parameter to below 0.1. The effective expo- 
sure after all screening was 89.6 ksec. Background data 



were parametrized using pcabackest v3 . 0, using the faint 
background model. Response matrices were created with 
pcarsp vlO.l, and nightly and total spectra were ex- 
tracted with saextrct v4.2d. Spectral fitting was per- 
formed with XSPEC vll.3.1, using PCA channels 5-32 
(corresponding approximately to 3-15 keV). To account 
for Galactic absorption, a column density of ATh = 
5.76 x 10 20 cm~ 2 w as used in the spectral fitting (cf., e.g., 
IWolter et all (|1998D . and also the PIMMS nH program). The 
influence of N-r is however marginal at this energy band. 

No flux variability was found within any of the 11 ob- 
serving nights. Between nightly averages, moderate changes 
were observed (the probability for a constant flux is 10 -5 , 
using the nightly averaged count rates). No hint for spectral 
variability was found. The light curve shown in Fig. [3] was 
derived by fixing the spectral model to the broken power 
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Tabs 


r*b 


Tdcabs 


r* 

x deabs 




0.23-4.0 TeV 


0.23-16.9 TeV 


0.23-4.0 TeV 


0.23-7.1 TeV 


All Data 
March 2005 
June 2004 


2 Q4+ ' 21 

3 1R+ - 48 

°- iD -0.61 


3-iol u ;i' 


L - 0± -Q.19 
1 4Q+0-19 

J -* a -o.ao 
1 7n+ - 47 

L - ' u -0.61 


1 74+U-ib 

x - 0.17 

1 R8+ 016 

2 1Q+ - 40 



Table 3. Photon indices from power-law fits to the VHE spectra of 1ES 1101-232. r a b s and r* bs are from fits to the 
measured spectra, Tdeabs and r*j eabs from fits to the deabsorbed spectra as described in Section 13.7.11 r a b s and Tdeabs 
correspond to the fits shown as solid lines in Fig.0J Photon indices T* bs and rj oabs are from fits including spectral 
bins above 4 TeV, which are compatible with zero flux, under the assumption that the power laws extend to these high 
energies. 



law derived for the total spectrum (see next paragraph), 
while leaving the normalisation as the only free parame- 
ter. Within statistical errors, the nightly fluxes derived are 
compatible with fluxes derived from single power-law fits 
with two free parameters (slope and normalisation). We 
conclude that 1ES 1101-232 was probably in a non-flaring 
state during this observation campaign. The simultaneous 
H.E.S.S. measurements are not sensitive to VHE 7-ray flux 
changes of similar variability amplitudes. 

As no spectral variability and only moderate flux 
changes were observed, a single averaged energy spectrum 
for the entire data set was derived. The spectrum be- 
tween 3 and 15keV is incompatible with a pure power 
law (Xrcd — 2.08 for 26 d.o.f., null-hypothesis probabil- 
ity 10 ), but a broken power law yields an acceptable 
fit with x^ ed = 1.18 (24 d.o.f., null- hypothesis probability 
25%). The total unabsorbed flux is i^-iokeV = (5-07 ± 
0.02 sta t) x 10~ 11 ergcm _2 s _1 . In Fig.[5j the unfolded spec- 
trum is shown. The spectrum is quite soft, with a pho- 
ton index of T = 2.49 ± 0.02 sta t below the break en- 
ergy inbreak = (7.9 ± 1.0) keV, and a slight softening to 
r = 2.78toii above -Ebreak- The spectrum shows that the 
peak in the SED was below ~ 3keV during the observa- 
tions. 

Possible systematic errors of the PCU spectrum were 
investigated by comparing recent archived spectra from 
Cassiopeia A with data from previous years, see the recom- 
mendations in the RXTE frequently asked question^. For 
the purpose of the analysis presented here, we concluded 
that systematic errors are of no concern. 

For the construction of the simultaneous SED discussed 
in Section [5j the X-ray spectrum derived from the entire 
RXTE data set was used, and similarly the H.E.S.S. spec- 
trum from the entire March 2005 data set. A restriction 
of both the H.E.S.S. and RXTE data to the strictly si- 
multaneous times appears overcautious, given the steady 
measured fluxes, and would have reduced the statistical ac- 
curacy, since only 13.3 (54%) of 24.9 hours RXTE observa- 
tions are overlapping with H.E.S.S. data. These 13.3 hours 
conversely represent 37% of the total March 2005 H.E.S.S. 
data set (35.6 hours of on-source observation time). We 
note that the March 2005 data set presents - in terms of si- 
multaneity and statistical accuracy - the best VHE/X-ray 
data set for 1ES 1101-232 so far. 
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Fig. 5. Unfolded X-ray spectra from 1ES 1101-232 in 
i?F(£')-representation. Points between 3-15 keV labeled 
Mar 5-16, 2005 are from RXTE, the line shows a bro- 
ken power-law fit to the data. Points labeled Jun 8 th , 
2004 between 0.2-10 keV are from XMM-Newton. The filled 
squares represent the reconstructed spectrum under the as- 
sumption of pure Galactic hydrogen absorption, and can be 
fit by a broken power law. For the spectrum shown with 
open squares, in addition to Galactic also absorption in the 
source was allowed, under the assumption of a pure power- 
law emission spectrum. 



4.1.2. Optical data 

During the March 2005 campaign, optical data on 
1ES 1101-232 were o btained using the ROTSE 3c telescope 
(jAkerlof et al.ll2003l ). which is located at the H.E.S.S. site. 
The ROTSE 3c telescope is a fast slewing robotic telescope 
with a 45 cm mirror and a wide field of view of 1.85°. To 
collect as much light as possible^, no optical filter is used; 
the ROTSE passband is 400 to 900 mn. 

During each of the 11 nights, typically 18 frames of 
60 seconds exposure time were obtained. After quality se- 
lection, 141 frames were used for analysis. The standard 



i.e., the probability that the assumed function fits the data 
\http: //heasarc. gsfc.nasa.gov/docs/xte/ftools/xtefaq. htrnl 



4 The main purpose of the ROTSE telescopes are optical 
afterglow observations of gamma ray bursts. 
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Satellite 


Observation ID 


Observation Dates 


Pointings 


Used Detectors And Live Time 


RXTE 


91123 


March 5-16, 2005 


48 


PCU 2: 89.60 ksec 


XMM-Newton 


205920601 


June 8, 2004 


1 


MOS2 (thin filter): 17.65 ksec 
PN (thin filter): 17.01 ksec 

OM: V 1.60, B 1.88, U 1.60, UVW1 4.00, UVM2 4.00 ksec 



Table 4. Observation log of the two X-ray satellite observations on 1ES 1101-232 used for this paper. 



ROTSE calibration chair0 delivered bias-corrected, flat- 
fielded images. To obtain object magnitudes, the stan- 
dard ROTSE procedure involves a vignetting correction 
and a relative photometry to a large number of compar- 
ison objects in th e field of view, using their USNO-A2 
(jMonet et al.lll998l ) R-band magnitudes as reference. The 
standard ROTSE photometry however failed to analyze the 
data from 1ES 1101-232 because of source confusion, there- 
fore a manual photometry was performed, as described in 
the following. 

Three temporally-stable comparison stars close to 
1ES 1101-232 were selecte d . Two of these stars have been 
identified by I Smith et alj (|l99lt ) as calibrators and have 
known photometric multi-color data, the third star was 
used for additional cross-checks. An aperture of R = 7.5" 
was used for photometry. The stability of our photome- 
try procedure was verified with several individual frames, 
by checking the correlation between measured fluxes and 
C-band magnitudes of the reference and various other 
field stars. C-band magnitudes (579-642 nm bandpass) 
were derived from th e CCD astrograph catalog UCAC2 
(jZacharias et al.ll2004h . 

To obtain a light curve and check for variability of the 
optical flux, for each frame an effective C-band magnitude 
of 1ES 1101-232 was derived, using the three comparison 
stars as calibrators. The optical light curve so derived ex- 
hibited only marginal flux variations, nightly flux averages 
showed changes below 0.1 mag. 

To derive an average optical flux of the BL Lac for use in 
the SED, an absolute flux calibration and subtraction of the 
host galaxy flux had to be performed. However, two facts 
caused this to be very difficult: the wide ROTSE bandpass, 
and the measured flux of m,Q — 16.03 mag which is close to 
the flux from the host galaxy (mj = 16.41 mag). Therefore, 
only an upper limit and a tentative lower limit of the optical 
flux could be derived. 

In order to subtract the flux from the host galaxy, it was 
verified that the s pectra of the two comparison stars used 
(Smi th et al.l!99i|) are similar t o template spectra of ellipti- 
cal hosts (jFukugita et aT]|l995l ) at the redshift of 1ES 1101- 
232. Applying the de Vaucouleurs profile with an effective 
radius of the host galaxy of R c = 4.1" and its tota l flux o f 
hir = 16.41 mag as measured by iFalomo fc Ulrichl |2000), 
we concluded that ~40% of the measured intensity comes 
from the host galaxy. If the BL Lac had also a similar spec- 
trum, the apparent magnitude of the BL Lac would then 
be mR = 16.4 mag. 

However, the wide bandpass of the ROTSE instrument 
causes this estimate to be quite susceptible to the actual BL 
Lac spectrum, which is expected to be much harder than 
the spectrum of the host galaxy. To estimate the magnitude 
of this effect, power-law spectra S v oc v~ a with a between 



5 http:/ / rotsel. physics. Isa. umich. edu, 



1.0 and -0.5 were simulated, under the assumption of a 
flat or symmetric response of the ROTSE detector between 
400 nm and 900 nm, and correction factors between 1.20 
and 1.44 to the R-band flux were derived. 

Magnitudes were finally corrected for Galactic extinc- 
tion, using a B-ban d extinction Ab = 0.254 mag (pro- 
vided by NED; from ISchleeel et all |1998D) and following 
the extinction laws of ICardelli et alj (fl989). Fluxes were 
derived from the ma gnitudes using the absolute calibra- 
tions bv lBessell (Il979l) . With this, an upper limit (assuming 
a = —0.5) of Sjx = 2.6 mjy was derived, and a tentative 
lower limit (using no correction factor) of Sr, = 1.0 mjy. 

4.2. Observations in June 2004 

The XMM-Newton observations on 1ES 1101-232 were con- 
ducted on June 8, 2004, as part of the pointings for XMM- 
Newton proposal 20592. A continuous pointing of 19.6 ksec 
was scheduled such that H.E.S.S. observations could be con- 
ducted simultaneously. Simultaneous ROTSE 3c observa- 
tions could not be performed because the telescope was 
under repair during that period. 

4.2.1. X-ray data 

The XMM-Newton data files were processed with 
xmmsas 6.5.0 following standard procedures as de- 
scribed in the SAS Handbook and calibration documents 
(jKirsch et al.1120061 ). where also systematic errors are dis- 
cussed. For our analysis, the most recent calibration files 
as available in January 2006 were used. Spectral and 
timing analysis was performed with XSPEC 11.3. Id and 
XR0N0S 5.21, respectively. 

The EPIC instruments during this pointing were set in 
timing (PN and MOS 1 cameras) and small window (MOS 2 
camera) mode, since one of the goals was to study spectral 
variability at the shortest possible timescales without pile- 
up problems, especially if the source were found in a very 
bright state. 

During this observation, however, the source was char- 
acterized by a constant flux on all timescales. Analysis of 
the light curves extracted in different energy bands, and 
the corresponding hardness ratios, showed no indication of 
spectral variability as well. Therefore the entire dataset was 
used to derive a time-integrated energy spectrum. In the 
following, since the MOS 1 data are in agreement with the 
other instruments but present a higher noise, the analysis 
was restricted to PN and MOS 2. 

Pile-up effects in the PN and MOS 2 data were checked 
with epatplot, and a mild pile-up was found for MOS 2, 
that can be removed considering single pixel events only 
(PATTERN=0), which were compared to spectra extracted 
with patterns 0-12 from different annular source regions 
dMolendi fc Sembavl l2003h . For the PN, since in timing 
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mode, PATTERN<4 was selected (single and double pixel 
events). For both instruments, quality FLAG=0 was used. 
The total livetime after screening for soft proton flares was 
17.005 ksec for PN and 16.653 ksec for MOS2. 

For MOS2, source counts were extracted from a cir- 
cle with radius 45" centered on the source centroid (ob- 
tained through eregionanalyse). The background was ex- 
tracted from the nearest source-free region on the periph- 
eral CCDs, with radius 90" . A check using different loca- 
tions on the peripheral CCDs, and the same source posi- 
tion on the blank-sky fields (available at the XMM-Newton 
Science Operating Center), showed no relevant differences, 
as was expected given the source rate (~ 8ctss _1 ) and 
spectrum. For the PN (which in timing mode has only one- 
dimensional imaging capabilities), photons were selected 
from rows 27 < RAWX < 47 (i.e., ±10 RAW pixels around 
the source strip), and 2 < RAWX < 18 for the background. 
To avoid the increased noise at low energies, the energy 
range for PN was restricted to 0.5-10 keV, while events 
down to 0.2 keV were accepted for MP S 2, as allowed by 
the latest calibration (|Kirsch et al.ll2006h . Spectra were re- 
binned with grppha to have at least 20 counts per channel. 
Response matrices and ancilliary files were produced with 
rmf gen and arf gen. 

The spectra were fitted with source models plus inter- 
stellar medium a bsorption, us i ng th e model phabs with 
abundances from I Wilms et al.l (|2000l ) . The absorbing col- 
umn density was fixed to the Galactic value Nn = 5.76 x 
10 20 cm -2 , but also a model with additional free absorption 
at the source (z=0.186) was tested. The host galaxy of 1ES 
1101-232 is not a normal elliptical ga laxy, the measured 
boxy isophotes (|Falortio fc Ulrichl [2000) may be indicative 
of extra absorbers, though in previous X-ray observations 
especially with BeppoSAX such possibility was statistically 
disfavoured (see Section |4~3")) . 

The PN and MOS 2 spectra were inspected separately, 
but finally the data from the two instruments were fitted 
together, with a free constant to allow for a different nor- 
malization between the two instruments (which remained 
within 3%). The results from the combined fit are reported 
in Tab.[5] The x 2 -values of the combined fits are relatively 
high, but to a large extent this is owing to residual cross- 
calibration uncertainties (though strongly improved with 
SAS 6.5) combined with large statistics. 

With Galactic absorption, a single power law does not 
provide an acceptable fit for the single detectors as well 
as for the combined fit (x 2 cd = 1.880 for 1547d.o.f.). 
A broken power-law model significantly improves all fits, 
with Xr C d = 1-195 (1545 d.o.f., null-hypothesis probability 
2 ■ 10 - *) for the combined fit. The hard (r : = 1.94 ± 0.01) 
and soft (r2 = 2.19 ± 0.01) spectral indices locate the X- 
ray peak of the SED at the break energy -Ebreak = 1.1 keV. 
Looking into the individual camera fits, the soft slopes are 
in agreement (r 2 , PN = 2.19±g;g| vs. r 2 , M oS2 = 2.21±g;g|), 
while break positions and hard slopes are slightly differ- 
ing (ri lPN = 1.98±g;S, Sbrcak.PN - 1.09±£° 2 keV, vs. 

ri,MOS2 = l-9li(X02> -E-break,MOS2 = 1.45^^2 keV); the 

respective fits are well acceptable for PN but slightly dis- 
favoured for MOS 2 (null-hypothesis probability: PN 5%, 
MOS 2 0.3%). Reasons for the discrepancies are the differ- 
ent fit ranges towards the soft end of the spectra, and the 
already mentioned calibration uncertainties. 



The combined broken power-law fit, considering Nn 
as a free parameter, yields a column density of N-r = 
5.8 ± 0.7 x 10 20 cm~ 2 which is in very good agreement with 
the Galactic value. To test whether additional absorption 
at the source is compatible with the data, a single power 
law model with Galactic absorption and free absorption 
at z=0.186 was investigated. The results are less favoured 
than the fits with pure Galactic absorption (null-hypothesis 
probability: PN: 0.7%, MOS 2: 10^ 15 , combined fit: 3-l(T n , 
extra absorption Nn ~ 3 x 10 20 cm -2 ), and provide more 
skewed residuals with an excess at lower energies less com- 
patible with the hypothesis of higher absorption. The dif- 
ference to the models with Galactic absorption is however 
not large, also considering the fact that the MOS 2 spec- 
trum still exhibits some unmodeled residua at low energies 
(0.4-0.5 keV), which could possibly be because of the men- 
tioned detector calibration uncertainties and/or imperfect 
modeling of the hydrogen absorption (even free abundances 
improve only slightly the residuals). 

To obtain the unfolded energy spectrum for the SED, 
the results of the combined fit were used. The residuals of 
the MOS 2 data to the respective MOS 2 model were mul- 
tiplied with the model derived for the PN. With this pre- 
scription, the absolute flux calibration from PN is trusted, 
while the full spectral energy range from MOS 2 can be 
used. Finally, the data were rebinned logarithmically in en- 
ergy. The unfolded spectrum is shown in Fig. [51 together 
with a similarly derived spectrum under the assumption of 
a pure power law and additional absorption at the source, 
for comparison. The model fits in Fig. [5] were derived from 
the unfolded spectra and are shown to indicate the size of 
the residuals. We note that the unfolded spectra do not 
strongly depend on the exact shape of the modeled spec- 
trum that is used in the unfolding procedure (using for 
example a pure power law instead of a broken power law 
yields a compatible spectrum within errors). 

To conclude, while small amounts of extra absorption 
at the source - which would soften the slope below 1 keV 
- cannot be excluded based on the XMM-Newton spectra 
alone, there is good evidence that the X-ray spectrum from 
1ES 1101-232 exhibited a peak in the SED at -1 keV during 
our observations, similar to earlier X-ray observations of 
this source. 

For the June 2004 SED discussed in the next section, the 
data of the entire H.E.S.S. June data set (i.e., ±3 days from 
the XMM-Newton pointing) were taken, in order to ob- 
tain a significant signal from the VHE data. Quantitatively, 
3.4 (66%) of the 5.1 hours (MOS 2) XMM-Newton obser- 
vations have simultaneous H.E.S.S. data. These 3.4 hours 
conversely represent 37% of the June 2004 H.E.S.S. data 
set (9.1 hours of on-source observation time). 

4.2.2. Optical data 

During the observations on 1ES 1101-232, the optical mon- 
itor (OM) onboard XMM-Newton took five exposures with 
five different filters, from the V to the UV band, with a win- 
dow in fast mode on our target. Details are given in Tab.0J 
As in the EPIC instruments, no variability was found in 
any OM exposure. Therefore, a photometric spectrum was 
extracted from all filters. 

OM data were processed with xmmsas 6.5.0. Only data 
from the imaging mode were used in the following. OM 
count rates were extracted using the point source analysis 
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procedure of the OM photometry program. An aperture of 
R = 6" was used for the source, background counts were 
extracted from an annulus with 7" < R < 12.5", for all 
filters. Count rates were converted into fluxes according to 
the prescriptions of the XMM-Newton watch-out pagefl 
The conversion factors for a white dwarf were used, as rec- 
ommended by the OM calibration scientist (Nora Loiseau, 
priv. comm.). 

The point spread function of the OM is considerably 
wider in the UV (UVW1 and UVM2 filters) than in the 
other filters, requiring specific aperture corrections in those 
bands. The default aperture for these filters could not be 
used because of bright sources in the corresponding back- 
ground annulus, therefore the same source and background 
apertures were used as for the other filters. The validity 
of the aperture correction was tested by selecting different 
source and background aperture sizes. The systematic error 
on fluxes derived with the UV filters was estimated to be 
-20%. 

In the V and B bands, a contribution from the 
host galaxy needs to be corrected. Using the same de 
Vaucouleurs profile as in Section 14.1.21 61% and 15% of 
the host galaxy flux were estimated to fall into the source 
and background apertures, respectively. These fractions of 
the host galaxy flux were subt racted, using its R - band mag- 
nitude tor = 16.41 mag from lFalomo fc Ulrichl (I2000T) and 
the ell iptical galaxy spectral template from iFukugita et al.l 
(|1995l) . 

Fluxes were finally corrected for Galactic extinction, 
again using As = 0.254 mag. The derived optical spectrum 
is shown in the SED in Fig.[SJ lower left panel. Error bars 
at the two UV flux points are from our systematic uncer- 
tainty estimate. We note that the U-filter flux is probably 
the most reliable flux estimate for the BL Lac, because of 
the mentioned uncertainties in the UV filters and because 
of the host galaxy influence at larger wavelengths. 

4.3. Previous X-ray observations of 1ES 1101-232 

1ES 1101-232 has already been observed in previous 
point ed X-ray obs ervations with RO SAT (jWolter et alj 
1998 ), BeppoSAX (IWolter et al.l l2000D and XMM-Newton 
( Perlman et al.ll2005f ). Spectra and fluxes were significantly 
different comparing individual observations, but the flux 
changes were not large (±25%). For example, the two spec- 
tra taken with B eppoSAX in 1997 an d 1998 (dubbed as high 
and low state in lWolter et alj ((2000)) only showed a signif- 
icant change above the peak at — 1.4 keV, with a total flux 
difference of 50%. Results from the different observations, 
including the two observations discussed in this paper, are 
summarized in Tab. [5] 

In all observations, pure intrinsic power laws were re- 
jected if only Galactic absorption values were assumed. 
All observations could be modeled with a pure power law, 
under the assumption of extra absorption of the order of 
Nh ~ 3 x 10 20 cm -2 . Nevertheless, in cases where the sta- 
tistical accuracy was sufficient (SAX I, SAX II, XMMII, cf. 
Tab. [5]), significantly better fits were derived, using only 
Galactic absorption and a curved intrinsic spectrum (usu- 
ally modeled by a broken power law). 

Assuming hence no extra absorption at the source, the 
low-energy peak in the SED of 1ES 1101-232 has in most 



observations been located in the — 0.5—3.5 keV range. From 
the difference of the position of the sp ectral break derived 



from the data set SAX I (see Tab.ED) bylWolter et"aH (2000) 
on the one hand and by iDonato et al.l ( 2005J) on the other 



hand, and also following the discussion in fPerlman et alj 
(120051) . we conclude that a broken power law is only an 
approximation of the real spectrum. Nevertheless, this does 
not affect the SED modeling presented in the next section, 
as the used unfolded XMM-Newton X-ray spectrum does 
not significantly change if the broken power law model is 
replaced by a curve d model, such as the one described in 
iFossati et all (|2000T) . 

5. Discussion 

5.1. Constraints on the extragalactic background light 

The detection of VHE 7-ray emission from 1ES 1101-232 
was used to constrain the density o f the EBL flux in the 
wavelength range of — 0.75-3.5 /j,m (|AHA06al and Section 
13. 8j) . The measured average photon index of T = 2.9 ei- 
ther precludes high EBL levels, or indicates a much harder 
instrinsic spectrum than seen or expected in other VHE 
blazars. For the purpose of this paper, the idea was adopted 
that the intrinsic VHE 7-ray spectrum of 1ES 1101-232 is 
restricted to TvHEjdeabs > 1.5. This limit is empirically de- 
rived in all blazars of lower redshift (where EBL corrections 
are less severe). It is also theoretically expected in present 
standard leptonic or hadronic scenarios for blazar emis- 
sion e mploying shock acceleration models (jMalkov &: Drurvl 
|2001[ ). because the hardest energy index obtained for accel- 
erated particles is p = 1.5 which results in T > 1.5 for 
all emissio n processes which can realistically be assumed 
()AHA06al L As shown in lAHA06al the limit of r V HE,dcabs > 
1.5 applied to the measured spectrum of 1ES 1101-232 re- 
sults in an EBL upper limit of vF v = 14nWm~ 2 sr _1 at 
1.5 ,um for typical EBL shapes. 

We note that harder photon spectra result if the low en- 
ergy end of the particle distribution is altered from the as- 
sumptions mentioned above. A variety of possibilities have 
been discussed in the literature. Relativistic Maxwellian 
(pileup) particle energy spectra, produced in a balance of 
stochastic acceleration and ra diative cooling fS chlickeiserl 
1985; iHenri fc Pelletien Il991f ). have previ ously been in- 
voked in the context of blazar spectra (jSauge fc Henril 
I2004T ) . Radiative cooling of an injection spectru m flatter 
than 7~ 2 could also r esult in a pileup spectrum fPinkau 
119801 : lKardashevl [l962'). Also bulk motion Comptonization 
of ambient low-frequency photons by a cold unshocked ul- 
trarelativistic jet with a very large Lorentz factor of the 
order of 10 6 ..10 7 could yield a narrow, line-like spectrum 
(lAharonian et al.ll2002a[). 

Katarzvhski et alj (|2006h us ed the total VHE 7- 



http://xmm.vilspa.esa.es/sas/new/watchout 



spectrum of 1ES 1101-232 from lAHA06al . together with 
archival X-ray data that were not taken simultaneously 
with the H.E.S.S. data, and showed that emission from a 
power-law type particle spectrum with a low energy cut- 
off at large Lorentz factors ~ 10 5 can produce an intrinsic 
VHE 7-spectrum as hard as oc v + 3. This would allow for 
an EBL level of vF v (l.§nm) ~ 20nWm- 2 sr- 1 and fit the 
VHE data. We note that our simultaneous data show no 
indication for such effects in the synchrotron branch. For 
the purpose of this paper, we refrain from discussing such 
spectra further. 
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Table 5. Spectral fit results of selected X-ray observations on 1ES 1101-232. (a-e): Values from spectral fits with low 
probability: (a):0.7%. (b) : <2% (c): <0.5%. (d):<l ( r 9 (e ):«lQ- 10 . t : Reanalysis of the data set "SAX I" presented 
originally by IWolter et al.l (|2000T ). tt: iDonato et all (|2005f) used Nn <ga ,i = 6.05 x 10 20 cm~ 2 and additional absorption at 
the source (2.4+I; 2 , x 10 20 cm~ 2 ); for simplicity, the sum of these two values is quoted here. Ti, inbreak, T 2 from a broken 
power-law fit, except * : Ti = T0.5-1.0kcV, T 2 = T4.0-10.0keV- The XMM II values are from the combined TN/MOS2 fits 
as described in the text, the flux normalisation is from TN only. The XTE row is detached to emphasi ze that the data are 
not sensitive to the br eak seen by all other in str uments in the ~ 0.5- 3.5 keV range. References: (1) IWolter et al.l fl 998: 
(2) IWolter et~aDl2000t (3) IDonato et~aTll2005l : (4) iPerlman et al.ll2005l . 



5.2. Spectral energy distribution of 1ES 1101-232 

For the construction of the SED, the deabsorbed data were 
taken, using what we consider the best available deabsorp- 
tion prescription. The optical and X-ray data were cor- 
rected for Galactic absorption, see Sections 14.1.21 ET2"TT1 and 
14.2.21 The VHE 7-ray spectra were deabsorbed, using two 
possible levels of the present day EBL photon field: the 
maximum value of vF v = 14nWm _2 sr _1 at 1.5 /im, and 
a minimum value of 10nWm _2 sr _1 co rresponding to the 
lower limit placed by Galaxy counts (jMadau fc Pozzettil 
2000). We note that the galaxy counts are presumably 
below the actual EBL density because of missed light 
(jTotani et al.l l200lh . To derive the optical depth for the 
VHE 7-r ays, the phenomenological EBL curve as used in 
AH A06al was applied, after scaling to match the given max- 
imum and minimum EBL densities and to take galaxy evo- 
lution effects into account, as explained in Section 1331 A 
redshift correction of the frequencies shown in the SED (to 
account for the difference of apparent and restframe wave- 
length) was not performed, as the correction would be dom- 
inated by the uncertainty of the emitting region's Doppler 
factor. 

In Fig.© the SEDs of 1ES 1101-232 are shown for the 
two periods for which broadband data together with the 
H.E.S.S. VHE measurements have been obtained. The up- 
per panel of Fig. [5] contains the average fluxes of the March 
2005 campaign (i.e., March 5-16, 2005), as derived from 
H.E.S.S RX TE, and ROTSE 3c data. As described in 
Section 14.1.11 the data have not been restricted to true si- 
multaneity (i.e., on minutes timescale). However, all data 
were taken during the same observing nights, with no sig- 
nificant (VHE) or only mild (X-ray, optical) variations be- 
tween days, and there is no sign of variability in these data 
on timescales shorter than a day. 

The SED shown in the lower panel of Fig. [5] has a lesser 
degree of simultaneity, because it contains the average VHE 
7-ray spectrum obtained from the H.E.S.S. June 2004 ob- 
servations (i.e., June 5-10, 2004), together with the X-ray 
and optical data as derived from the XMM-Newton obser- 
vation performed on June 8 (see Section 14.2. 1[) . 



5.3. Model-independent considerations 

Independent of specific emission scenarios, the broadband 
data from 1ES 1101-232 presented here show some remark- 
able features: 

VHE peak: The hard intrinsic VHE 7-ray spectrum 
(Tdoabs S 2) requires that the VHE 7-ray flux peak is lo- 
cated at > 3 TeV and that the measured VHE photons do 
come from below that peak, unless the lowest possible EBL 
level is assumed and errors on the VHE 7-ray spectrum are 
exploited to the limits. The VHE 7-ray spectrum does not 
show hints of curvature. A steepening towards higher en- 
ergies would indicate the position of the SED high-energy 
peak at ~ 3 TeV (either due to a break in the particle dis- 
tribution, or because of Klein Nishina effects) , but the data 
are insufficient to claim such a steepening. 

Comparison 2004/2005: The source was in a different 
state in March 2005 compared with the June 2004 period, 
judging from the X-ray data. While the X-ray flux and spec- 
trum in June 2004 were comparable to previous X-ray mea- 
surements (see Tabj5]), the X-ray spectrum in March 2005 
as measured with RXTE shows a higher flux and a fairly 
soft spectrum, softer than any previously measured X-ray 
spectrum from this source. The synchrotron peak (usually 
at ~ 0.5-3.5 keV) cannot have shifted to higher energies 
during the March 2005 high state, the XTE spectrum pre- 
cludes a peak above ~ 3keV. Despite the different X-ray 
flux states in March 2005 and June 2004, there is no evi- 
dence for a change of the VHE 7-ray spectrum between the 
two periods; nevertheless, statistical errors from the June 
2004 data sets could allow for a factor of up to about two 
flux difference between the two periods. 

X-ray variability: The X-ray light curve in March 2005 
only exhibits mild variations. Under the assumption of 7~ 2 - 
type injection spectra, it seems unlikely that the soft RXTE 
spectrum from that period is governed by emission from a 
cooled particle spectrum. The spectrum could reflect the 
cut-off of the acceleration spectrum, but depending on the 
assu med scenario, other op tions like softer injection spectra 
(e.g.. iGhisellini et al.ll2002f ) are also possible. 

Leptonic emission spectrum: One can compare X-ray 
and VHE 7-ray spectral indices under the assumption that 
electrons from the same energy band with iV 7 ex 7 _p are 
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Fig. 6. Spectral energy distribution of 1ES 1101-232. Upper panel: Data from March 5-16, 2005. X-ray data are from 
RXTE. In the optical band, an upper limit (filled triangle) and a tentative lower limit (open triangle) from ROTSE 3c 
data are shown, see Section 14.1.21 for details. In the VHE band, the measured H.E.S.S. spectrum (red, open symbols) 
and the deabsorbed spectrum using a maximum EBL level of 14nWm _2 sr -1 at 1.5 (im (see text) are shown; for better 
visibility, all open symbols were slightly shifted to the left, to 90% of the respective original frequency. The thick dashed 
line is a power-law fit to the deabsorbed data as plotted, while the thin dashed line indicates the effect if the EBL is 
lowered to the minimum level of 10nWm~ 2 sr -1 . The latter value corresponds to the minimum level of the EBL as derived 
from known resolved galaxy counts. Thick and thin solid curves denote results from a single zone SSC model. The thick 
curves represent a model that was optimized to fit the H.E.S.S. and X-ray data, while the thin lines denote a model with 
an electron distribution A r . oc j~ 2 below the break. Lower panel: Data from June 2004. X-ray and optical data were 
derived from an XMM-Newton pointing on June 8, 2004. In the VHE band, H.E.S.S. data taken between June 5-10, 
2004, are shown, using the same procedure as described for the upper panel. Solid curves denote results from a single 
zone emission model, also following the same prescription as for the upper panel. 



emitting X-rays via synchrotron and VHE photons via IC 
processes. The rather flat X-ray spectrum over nearly two 
decades in energy, as seen with XMM-Newton in 2004, con- 
strains p to 2.9 < p < 3.4, and therefore requires that 
the VHE 7-ray emission is in the Thompson regime (1.9 < 
Tic < 2.2) to be still in agreement with the measured VHE 
7-ray spectrum (rj un 2oo4 = 1-70 ± 0.47). Assuming that 
the VHE 7-ray spectrum was constant throughout the ob- 
servation periods (r to tai = 1.51 ±0.17) would require either 
to drop the initial assumption about the common spectral 
range of the emitting electrons, or to lower the EBL to the 
lowest possible value of vF v (1.5 /xm) = 10nWm _2 sr _1 . For 
the March 2005 data set, a common energy band of the 
electron spectrum cannot account for the X-ray and VHE 
7-ray emission. 

Intrinsic source luminosity: Depending on the assumed 
emission mechanism and EBL density, estimates of the in- 
trinsic luminosity Lvhe of the VHE 7-ray emitting region 



in 1ES 1101- 2 32 can reach unusually large values. Following 
Siko ra et al.l (|1997[ ). the intrinsic luminosity would be 
Lvhe > 10 39 ergs _1 under the assumption of SSC VHE 7- 
ray emission from a moving blob with Doppler factor (5^30. 
Since however the emission seems constant on timescales of 
days or maybe even years, a stationary jet model could also 
be invoked, leadi ng to estimates of Lvhe > 10 42 ergs~ 1 . 
Adopting like in iKatarzvnski et aTl (l2006f) the "best fit" 
EBL model from iKneiske et all (120041) . with an EBL den- 



sity approaching uF v (1.5 /im) ~ 20nWm 2 sr , would in- 



crease the intrinsic luminosity estimates to > 10 40 ergs" 
for the moving blob and even > 10 43 ergs -1 for the station 
ary scenario. 



5.4. One-zone leptonic emission model 

Neither flux correlations between the X-ray and VHE bands 
nor fast variability (i.e., on sub-day or even sub- hour 
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2005 I 


2005 II 


2004 I 


2004 II 


Pi 


2 


1.5 


2 


1.7 


P2 


4.3 


7.0 


3.6 


3.6 


7b 


1.8 ■ 10 s 


3.8 • 10 5 


1.8 ■ 10 5 


1.8 • 10 5 


A"[cm- 3 ] 


3.5 • 10 2 


15 


9.0 ■ 10 2 


40 


R[cm] 


2.8 ■ 10 16 


5.7- 10 15 


1.7 • 10 16 


1.15 ■ 10 16 



8 25 
B[G] 0.1 

7min 10 3 
Tmax 10° 



Table 6. Physical parameters of the SSC one-zone model 
spectra. 2005 I and 2004 I, i.e. the scenarios with p\ = 2, 
correspond to the models shown as thin solid lines in Fig. [5] 
2005 II and 2004 II correspond to the models shown as thick 
solid lines. Doppler factor, magnetic field, 7 m i n , and 7 max 
are in all four cases identical (Hq — 75 kms _1 Mpc _1 ). 



timescale) could be established with the present data set. 
Therefore, there is no clear indication that would favour 
one-zone emissio n models, either leptonic or hadronic. Past 
experience (e.g., lAharonian et al.l l2005af ) has shown that 
both hadronic and leptonic scenarios have enough free pa- 
rameters to generally be able to explain the broad-band 
emission from blazars, if no further arguments from vari- 
ability can be invoked. Nevertheless, the reconstruced hard 
VHE 7-ray spectrum from 1ES 1101-232 challenges one- 
zone emission models. 

A simple leptonic emission model was used, a time- 
indepe ndent SSC model as described bv lKatarzvhski et al.l 
(|200lh . in an attempt to describe the measured VHE, X-ray, 
and optical data. The description of a one-zone homoge- 
neous, spherical emitting region R and a homogeneous mag- 
netic field B was adopted, which propagates with Doppler 
factor 5 towards the observer. The high-energy electron 
distribution is modeled with a broken power law between 
Lorentz factors 7 m in and 7 max with a break at 7b, and a 
density normalisation K. The two epochs have been con- 
sidered independently. In the March 2005 data set, the most 
stringent constraint on the SED comes from the hard 7-ray 
spectrum extending towards 3 TeV. In the June 2004 case, 
the 7-ray part of the spectrum is easier to reproduce be- 
cause of the smaller number of events and the slightly softer 
spectrum, while the X-rays are more constraining with a 
rather flat spectrum over a large frequency range. 

The results from two parameter sets for each SED are 
shown in Fig.[5]and Tab. [51 In case I, p\ (the energy index of 
the particles between 7 m i n and 7b) was set to 2, as expected 
from an uncooled, shock-accelerated particle distribution. 
With 7b, 2005 = 1-8 • 10 5 and p 2 ,2005 = 4.3 (p 2 being the 
energy index between 7b and 7 max ) for the March 2005 
SED, and 7b, 2004 = 1-8 ■ 10 5 , ^2,2004 = 3.6 for the June 
2004 data, good fits were obtained to the X-ray and optical 
data, respectively. As expected, the hard VHE spectra are 
not well reproduced when using p = 2. The 2004 VHE 
data are still satisfactorily matched, but the very hard 2005 
VHE 7-ray spectrum (reconstructed with the EBL density 
14 nWm _2 sr _1 ), having also smaller statistical errors, is not 
well reproduced by the model. A simple \ 2 test only yields a 
just over 2 a compatibility between the data and the model. 
A lower EBL level improves the agreement with the data. 



In case II, p\ was chosen to be of the order of 1.5, which 
can be expected for instance from particle acceleration at 
strong shocks in a relativistic gas. Better fits are then ob- 
tained for the 7-ray spectra at both epochs as illustrated 
in Fig. [SI especially for 2005 where the \ 2 test now yields a 
1 a compatibility. In this scenario, the optical flux cannot 
be described within the SSC one-zone model, and has to 
be attributed to another component. This additional low 
frequency emission could come for example from some ex- 
tended jet emission. Such an additional component is in 
any case needed to explain the emission at radio frequen- 
cies, which was m easured with the Nancay radio telescope 
at 2.685 GHz (see lAharonian et al.l (|2005al ) for a descrip- 
tion of the instrument and data processing). The flux lev- 
els obtained in June 2004 (0.11 ± 0.02 Jy, not simultaneous 
to the other data presented in this paper) and in March 
2005 (0.08 ± 0.01 Jy, simultaneous to the data presented 
i n this paper) were comparable to previ ous measurements 
(|Griffith et al . I H99l [Douglas et al]ll996h . 

To conclude, SSC one-zone models are globally able to 
reproduce the SED of 1ES 1101-232 from the X-rays to 
VHE 7-rays for the data set analyzed here. For the 2005 
data set, an inclusion of the optical lower limit, obtained 
through the ROTSE 3c detector, is however statistically 
disfavoured. Moreover we should stress that with 1ES 1101- 
232 the limit of the capabili ties of SS C one-zone models is 
reached, as also discussed in lAHA06al It is very difficult to 
get good fits for the shape of the VHE 7-ray tail of the ob- 
served spectra, as long as one keeps usual assumptions for 
particle acceleration mechanisms and does not assume an 
EBL level as low as 10nWm _2 sr" 1 . The generated spectra 
deviate, especially for the March 2005 data, from the hard 
spectra obtained with H.E.S.S. Smaller statistical errors on 
the VHE 7-ray spectrum or an extension to higher energies 
(or both), further constraints on the size of the emitting 
zone from variability data, or a slight increase of the ab- 
sorption by extragalactic background above the value of 
14nWm _2 sr _1 , may reach the limit of one-zone SSC mod- 
els. More complex (e.g., two-zone) scenarios may therefore 
be required. 

6. Conclusion 

Using the H.E.S.S. Cherenkov telescope system, VHE 7-ray 
emission from 1ES 1101-232 was discovered. 1ES 1101-232 
is currently the farthest object (z=0.186) with confirmed 
redshift detected in VHE 7-rays. The VHE 7-ray spectrum 
favours a low level of the extragalactic background light at 
a few microns, close to the lower limit placed by galaxy 
counts. 

The VHE data from 1ES 1101-232 were taken together 
with measurements in the X-ray and optical bands. The 
best SED from truly simultaneous observations could be 
constructed from a multifrcqucncy campaign performed in 
March 5-16, 2005. The data revealed that the source was 
brighter in X-rays during this campaign than in any other 
previous X-ray observation, but did not show strong flaring 
activity in either the X-ray or the VHE band. The H.E.S.S. 
data set is compatible with constant VHE 7-ray emission 
throughout all observation periods. 

Using an EBL density of vF v — 14nWm _2 sr _1 at 
1.5 /mi to deabsorb the VHE 7-ray spectrum, the H.E.S.S. 
data result in a very hard intrinsic spectrum of 1ES 1101- 
232, with a peak in the VHE power output above 3 TeV. 
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The spectrum is harder (T ~ 1.5) than in parametrisations 
using SSC models. An EBL level below 14nWm~ 2 sr _1 
would soften the 7-ray spectrum nearly to T ~ 2.0, which 
would ease the application of SSC scenarios but at the 
same time challenge EBL models. On the other hand, an 
EBL density above 14nWm _2 sr _1 would result in an even 
harder 7-ray spectrum, therefore such high EBL levels can 
be excluded under the assumption that 1ES 1101- 232 is sim- 
ilar to all VHE blazars known so far (AHA06a). A 7-ray 
spectrum in 1ES 1101-232 harder than T ~ 1.5 would indi- 
cate previously unseen blazar physics, but the broadband 
data presented in this paper currently do not support such 
a conclusion. 

In conclusion we find 1ES 1101-232 to be at an extreme 
end of blazar characteristics. This is the first time a de- 
tected blazar spectrum peaks above the VHE band. The 
data challenge current blazar models in the high-energy end 
of the electromagnetic spectrum. 
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